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Abstract 

We use distributed computing tools to provide a new perspective on the behavior of 
cooperative biological ensembles. We introduce the Ants Nearby Treasure Search (ANTS) 
problem, a generalization of the classical cow-path problem [101 1201 1411 W2\ , which is relevant 
for collective foraging in animal groups. In the ANTS problem, k identical (probabilistic) 
agents, initially placed at some central location, collectively search for a treasure in the 
two-dimensional plane. The treasure is placed at a target location by an adversary and the 
goal is to find it as fast as possible as a function of both k and D, where D is the distance 
between the central location and the target. This is biologically motivated by cooperative, 
central place foraging, such as performed by ants around their nest. In this type of search 
there is a strong preference to locate nearby food sources before those that are further away. 
We focus on trying to find what can be achieved if communication is limited or altogether 
absent. Indeed, to avoid overlaps agents must be highly dispersed making communication 
difficult. Furthermore, if the agents do not commence the search in synchrony, then even 
initial communication is problematic. This holds, in particular, with respect to the question 
of whether the agents can communicate and conclude their total number, k. It turns out 
that the knowledge of k by the individual agents is crucial for performance. Indeed, it is a 
straightforward observation that the time required for finding the treasure is Q,(D + D 2 /k), 
and we show in this paper that this bound can be matched if the agents have knowledge of k 
up to some constant approximation. 

We present a tight bound for the competitive penalty that must be paid, in the running 
time, if the agents have no information about k. Specifically, this bound is slightly more 
than logarithmic in the number of agents. In addition, we give a lower bound for the setting 
in which the agents are given some estimation of k. Informally, our results imply that the 
agents can potentially perform well without any knowledge of their total number k, however, 
to further improve, they must use non trivial information regarding k. Finally, we propose a 
uniform algorithm that is both efficient and extremely simple, suggesting its relevance for 
actual biological scenarios. 

1 Introduction 

1.1 Background and Motivation 

The universality of search behaviors is reflected in multitudes of studies in different fields including 
control systems, distributed computing and biology. We use tools from distributed computing to 
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study a biologically inspired scenario in which a group of agents, initially located at one central 
location, cooperatively search for treasures in the plane. The goal of the search is to locate nearby 
treasures as fast as possible and at a rate that scales well with the number of participating agents. 

A variety of animal species search for food around a central location that serves the search's 
initial point, final destination or both [17] ■ This central location could be a food storage area, 
a nest where offspring are reared or simply a sheltered or familiar environment. Central place 
foraging holds a strong preference to locating nearby food sources before those that are further 
away. Possible reasons for that are, for example: (1) decreasing predation risk |44:, (2) increasing 
the rate of food collection once a large quantity of food is found [57], (3) holding a territory 
without the need to reclaim it [35] HU H5J , and (4) the ease of navigating back after collecting the 
food using familiar landmarks [T5] . 

Searching in groups can increase foraging efficiency 38, p. 732]. In some extreme cases, food 
is so scarce that group searching is believed to be required for survival |18( I39j . Proximity of the 
food source to the central location is again important in this case. For example, in the case of 
recruitment, a nearby food source would be beneficial not only to the individual that located the 
source but also increase the subsequent retrieval rate for many other collaborators [50]. Foraging 
in groups can also facilitate the defense of larger territories [3S]. Eusocial insects (e.g., bees 
and ants) engage in highly cooperative foraging, which can be expected as these insects reduce 
competition between individuals to a minimum and share any food that is found. Social insects 
often live in a single nest or hive, which naturally makes their foraging patterns central. 

Little is known about the communication between foragers, but it is believed that in some 
scenarios communication may become impractical |36j . This holds, for example, if the foragers 
start the search at different times and remain far apart (which may be necessary to avoid 
unnecessary overlaps). Hence, the question of how efficient the search can be if the communication 
is limited, or altogether absent, is of great importance. 

In this paper, we theoretically address general questions of collective searches in the particular 
natural setting described above. More precisely, we introduce the Ants Nearby Treasure Search 
(ANTS) problem, in which k identical (probabilistic) agents, initially placed at some central 
location, collectively search for a treasure in the two-dimensional plane. The treasure is placed 
by an adversary at some target location at distance D from the central location, where D is 
unknown to the agents. The goal of the agents is to find the treasure as fast as possible, where 
the time complexity is evaluated as a function of both k and D. 

In the context of search algorithms, evaluating the time as a function of D was first introduced 
in a classical paper by Baeza- Yates et al. [TO], who studied the cow-path problem (studied also 
in [2"0 ] I1T | I4"2"]). The ANTS problem generalizes the cow-path problem, as it considers multiple 
identical agents instead of a single agent (a cow in their terminology). Indeed, in this distributed 
setting, we are concerned with the speed-up measure (see also, [5] M H7J HH]), which aims to 
capture the impact of using k searchers in comparison to using a single one. Note that the 
objectives of quickly finding nearby treasures and having significant speed-up may be at conflict. 
That is, in order to ensure that nearby treasures are quickly found, a large enough fraction of the 
search force must be deployed near the central location. In turn, this crowding can potentially 
lead to overlapping searches that decrease individual efficiency. 

It is a rather straightforward observation that the time required for finding the treasure is 
fl(D + D 2 /k). Our focus is on the question of how agents can approach this bound if their 
communication is limited or even completely absent. In particular, as information of the size of 
the foraging group may not be available to the individual searchers, we concentrate our attention 
on the question of how important it is for agents to know (or estimate) their total number. As 
we later show, the lack of such knowledge may have a non-negligible impact on the performance. 
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1.2 Our Results 



We introduce and investigate the ANTS problem, a generalization of the cow-path problem. First, 
we show that if the agents have a constant approximation of their total number k, then there 
exists a rather simple search algorithm whose expected running time is 0(D + D 2 /k), making 
it 0(l)-competitive. We then turn our attention to uniform searching algorithms, in which the 
agents are not assumed to have any information regarding k. We completely characterize the 
speed-up penalty that must be paid when using uniform algorithms. Specifically, we show that 
for a function / such that Yl'jLi V/(i) converges, there exists a uniform search algorithm that is 
0(/(log fc))-competitive. On the other hand, we show that if Yl'jLi V/C?) diverges, then there is 
no uniform search algorithm that is 0(/(log fc))-competitive. In particular, this implies that, for 
every constant e > 0, there exists a uniform search algorithm that is 0(log 1+e fc)-competitive, 
but there is no uniform search algorithm that is 0(log /c)-competitive. Hence, the penalty for 
using uniform algorithms is slightly more than logarithmic in the number of agents. 

In addition, we give a lower bound for the intermediate setting in which the agents are given 
some estimation of k. As a special case, this lower bound implies that for any constant e > 0, 
if each agent is given a (one-sided) fc e -approximation of k, then the competitiveness is O(logfc). 
Informally, our results imply that the agents can potentially perform well without any knowledge 
of their total number k, however, to further improve they must use non trivial information 
regarding k. Finally, we propose a uniform search algorithm that is concurrently efficient and 
extremely simple and, as such, it may imply some relevance for actual biological scenarios. 

1.3 Related Work 

Our work falls within the scope of natural algorithms, a recent attempt to investigate biological 
phenomena from an algorithmic perspective [TJ 1131 115j . 

Collective search is a classical problem that has been extensively studied in different fields 
of science. Group living and food sources that have to be actively sought after make collective 
foraging a widespread biological phenomenon. Social foraging theory |33j makes use of economic 
and game theory to optimize food exploitation as a function of the group size and the degree of 
cooperativity between agents in different environmental settings. Social foraging theory has been 
extensively compared to experimental data (see, e.g., [5] EH]) but does not typically account for 
the spatial characteristics of resource abundance. Central place foraging theory [27] assumes a 
situation in which food is collected from a patchy resource and is returned to a particular location, 
such as a nest. This theory is used to calculate optimal durations for exploiting food patches 
at different distances from the central location and has also been tested against experimental 
observations |37] . Collective foraging around a central location is particularly interesting in 
the case of social insects where large groups forage cooperatively with, practically, no competition 
between individuals. Harkness and Mardouras |36j have conducted a joint experimental and 
modeling research into the collective search behavior of non-communicating desert ants. Modeling 
the ants' trajectories using biased random walks, they reproduce some of the experimental findings 
and demonstrate significant speed-up with group size. In bold contrast to these random walks, 
Reynolds [51] argues that Levy flights with a power law that approaches unity is the optimal 
search strategy for cooperative foragers as traveling in straight lines tends to decrease overlaps 
between searchers. 

From an engineering perspective, the distributed cooperation of a team of autonomous agents 
(often referred to as robots or UAVs: Unmanned Aerial Vehicles) is a problem that has been 
extensively studied. These models extend single agent searches in which an agent with limited 
sensing abilities attempts to locate one or several mobile or immobile targets |49j . The memory 
and computational capacities of the agent are typically large and many algorithms rely on the 
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construction of cognitive maps of the search area that includes current estimates that the target 
resides in each point |64j . The agent then plans an optimal path within this map with the intent, 
for example, of optimizing the rate of uncertainty decrease [40 . Cooperative searches typically 
include communication between the agents that can be transmitted up to a given distance, or even 
without any restriction. Models have been suggested where agents can communicate by altering 
the environment to which other agent then react |61j . Cooperation without communication has 
also been explored to some extent [7] but the analysis puts no emphasis on the speed-up of the 
search process. In addition, to the best of our knowledge, no works exist in this context that put 
emphasis on finding nearby targets faster than faraway one. Similar problems studied in this 
context are pattern formation [THl \±M , rendezvous [21I2B], and flocking [35]. It is important 
to stress that in all those engineering works, the issue of whether the robots know their total 
number is typically not addressed, as obtaining such information does not seem to be problematic. 
Furthermore, in many works, robots are not identical and have unique identities. 

In the theory of computer science, the exploration of graphs using mobile agents is a central 
question. Most of the research for graph exploration is concerned with the case of a single 
deterministic agent exploring a finite graph (typically, with some restrictions on the resources 
of the agent and/or on the graph structure). For example, in [3, 12K2IJ122] t ne agent explores 
strongly connected directed finite graphs, and in [5J [22] 12H I2H Ell HFJ E3] the agent explores 
undirected finite graphs. When it comes to probabilistic searching, the random walk is a natural 
candidate, as it is extremely simple, uses no memory, and trivially self-stabilizes. Unfortunately, 
however, the random walk turns out to be inefficient in a two-dimensional infinite grid. Specifically, 
in this case, the expected hitting time is infinite, even if the treasure is nearby. 

Evaluating the time to find the treasure as a function of D, the initial distance to the 
treasure, was studied in the context of the cow-path problem. One of the first papers that 
studied the cow-path problem is a paper by Baeza- Yates et al. [TO], in which the competitive 
ratio for deterministically finding a point on the real line is shown to be nine. Considering the 
two-dimensional case, Baeza- Yates et al. proved that the spiral search algorithm is optimal up to 
lower order terms. Randomized algorithms for the problem were studied by Kao et al. [5T] , for the 
infinite star topology. Karp et al. [12] studied an early variant of the cow-path problem on a binary 
tree. Recently, Demaine et al. |20j considered the cow-path problem with a double component 
price: the first is distance and the second is turn cost. Lopez-Ortiz and Sweet [55] extended 
the cow-path problem by considering k agents. However, in contrast to the ANTS problem, the 
agents they consider are not identical, and the goal is achieved by (centrally) designing a different 
specific path for each of the k agents. 

In general, the more complex setting of using multiple identical agents has received much less 
attention. Exploration by deterministic multiple agents was studied in, e.g., [8] [23 120] ■ To 
obtain better results when using several identical deterministic agents, one must assume that the 
agents are either centrally coordinated or that they have some means of communication (cither 
explicitly, or implicitly, by being able to detect the presence of nearby agents). When it comes to 
probabilistic agents, analyzing the speed-up measure for /c-random walkers has recently gained 
attention. In a series of papers, initiated by Alon et al. [3], a speed-up of f2(fc) is established for 
various finite graph families, including, in particular, expanders and random graphs [41 1171 [55] . 
While some graph families enjoy linear speed-up, for many graph classes, to obtain linear speed-up, 
k has to be quite small. In particular, this is true for the two-dimensional n-node grid, where 
a linear speed up is obtained when k < 0(log 1_e n). On the other hand, the cover time of 
2-dimensional n-node grid is always f2(n/log k), regardless of k. Hence, when k is polynomial in n, 
the speed up is only logarithmic in k. The situation with infinite grids is even worse. Specifically, 
though the /c-random walkers would find the treasure with probability one, the expected time to 
find the treasure becomes infinite. 
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The question of how important it is for individual processors to know their total number has 
recently been addressed in the context of locality. Generally speaking, it has been observed that 
for several classical local computation tasks, knowing the number of processors is not essential [35] . 
On the other hand, in the context of local decision, some evidence exist that such knowledge may 
be crucial for non-deterministic distributed decision [3D]. 

2 Preliminaries 

We consider the ANTS problem, where k mobile agents (robots) are searching for a treasure on the 
two-dimensional plane. Each agent has a field of view bounded by a small positive constant, hence, 
for simplicity, we can assume that the agents are actually walking on the integer two-dimensional 
infinite grid G := Z 2 . All k agents starts the search from a central node s of G, called the source. 
An adversary locates the treasure at some node r of G, referred to as the target node; the agents 
have no a priori information about the location of t. The goal of the agents it to find the treasure: 
this task is accomplished once at least one of the agents visits the node r. 

The agents are probabilistic machines that can move on the grid, but cannot communicate 
between themselves. All k agents are identical (execute the same protocol). An agent can traverse 
an edge of the grid in both directions. We do not restrict the internal storage and computational 
power of agents, nevertheless, we note that all our upper bounds use simple procedures that can 
be implemented using relatively short memory. For example, with respect to navigation, our 
constructions only assume the ability to perform four basic procedures, namely: (1) choose a 
direction uniformly at random, (2) walk in a "straight" line to a prescribed distance, (3) perform 
a spiral search around a nod^J and (4) return to the source node. On the other hand, for our 
lower bounds to hold, we do not require any restriction on the navigation capabilities. 

Regarding the time complexity, we assume that traversing an edge takes one unit of time. 
Furthermore, for the ease of presentation, we assume that the agents are synchronous, that is, 
each edge traversal costs precisely one unit of time (and all internal computations are performed 
in zero time). Indeed, this assumption can easily be removed if we measure the time according to 
the slowest edge-traversal. We also assume that all agents start the search simultaneously at the 
same time, denoted by t . This assumption can also be easily removed by starting to count the 
time after the last agent initiates the search. We measure the cost of an algorithm by its expected 
running time, that is, the expected time (from time to) until at least one of the agents finds the 
treasure. We denote the expected running time of algorithm A by Ta(D, k). 

The distance d(u, v) between two nodes u and v of G is simply the hop distance between 
them, i.e., the number of edges on a shortest path connecting u and v in G. The distance 
between a node u and the source node s is simply written d(u), i.e. d(u) := d(s,u). Let D 
be the distance between the source node s and the target node r, i.e., D := d(r). For a node 
u of G and an integer r, let B(u, r) be the ball of radius r centered at the node u, formally, 
B(u,r) :— {v G G : d(u,v) < r}. For short, we denote the ball of radius r around the source 
node s by -B(r), that is, B(r) := B(s,r). 

Note that if an agent knows D, then it can potentially find the treasure in time O(D), by 
going to any node at distance D, and then performing a circle around the source of radius D 
(assuming, of course, that its navigation abilities enable it to perform such a circle). On the other 

x Thc spiral search around a node v is a particular deterministic local search algorithm (see, e.g., |10|) that 
starts at v and enables the agent to visit all nodes at distance Q(y/x) from v by traversing x edges, for every 
integer x. For our purposes, since we are concerned only with asymptotic results, we can replace this atomic 
navigation procedure with any procedure that guarantees such a property. For simplicity, in the remainder of the 
paper, we assume that for every integer x, the spiral search of length x starting at a node v visits all nodes at 
distance at most \Jxj2 from v 
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hand, without knowledge about D, an agent can find the treasure in time 0(D) by performing a 
spiral search around the source. When considering k agents, it is easy to sec[^]that the expected 
running time is f2(D + D 2 /k), even if the total number of agents k is known to all agents, and 
even if we relax the model and allow agents to freely communicate between each other. It follows 
from Theorem [XT] that if k is known to agents then there exists a search algorithm whose expected 
running time is asymptotically optimal, namely, 0(D + D 2 /k). We evaluate the performance of 
an algorithm that does not assume the precise knowledge of k with respect to this aforementioned 
optimal time. Formally, let 4>(k) be a function of k. An algorithm A is <p(k)- competitive if 

T A (D,k) <cl>{k)-{D + D 2 /k), 

for every integers k and D. We shall be particularly interested in the performances of uniform 
algorithms: these are algorithms in which no information regarding k is available to agents. (The 
term "uniform" is chosen to stress that agents execute the same algorithm regardless of their 
number, see, e.g., [33].) 

3 Upper Bounds 

3.1 Optimal Running Time with Knowledge on k 

Our first theorem asserts that agents can achieve an asymptotically optimal running time if they 
know the precise value of k. As a corollary, it follows that, in fact, to obtain such a bound, it is 
sufficient to assume that agents only know a constant approximation of k. 

Theorem 3.1. Assume that the agents know their total number, k. Then, there exists a (non- 
uniform) search algorithm running in expected time 0(D + D 2 /k). 

Proof. For i e N, set Bi :— \u : d(u) < 2 1 }. Consider the following algorithm (see also Figure [l] 
for an illustration). 

Fix a positive integer I and consider the time Ti until each agent completed £ phases i with 
* _^ log D. Each time an agent performs phase i, the agent finds the treasure if the chosen node u 
belongs to the ball B(t, \ftl/2) around r, the node holding the treasure. Note that at least some 
constant fraction of the ball B(t, ^/ti/2) is contained in Bi, since i > log-D. The probability of 
choosing a node u in that fraction is thus VL(\B(t, y/Ti/2)\/\Bi\), which is at least (3/k for some 
positive constant /3. Thus, the probability that by time Te none of the k agents finds the treasure 
(while executing their respective I phases i) is at most (1 — f3/k) ke , which is at most 7 for some 
constant 7 greater than 1. 

For i 6 N, let ip(i) be the time required to execute a phase i. Note that ip(i) = 0(2 l + 2 2% /k). 
Hence, the time until all agents complete stage j for the first time is 

3 3 
53 ^(i) = 0{V +J2^ 2l /k) = 0(V + 2 2 Vfc). 

i=l i=l 

Now fix s := [log D~\ . It follows that for any integer £, all agents complete their respective stages 
s + I by time T(£) := 0(2 s+e + 2 2 ^ s+ ^ /k). Observe that by this time, all agents have completed 

2 To see why, consider a search algorithm A whose expected running time is T. Clearly, T > D, because it takes 
time D to merely reach the treasure. Assume, towards a contradiction, that T < D 2 /4k. In any execution of A, 
by time 2T, the k agents can visit a total of at most 2Tk < D 2 /2 nodes. Hence, by time 2T, more than half of the 
nodes in Bu '■= {u : 1 < d(u) < D} were not visited. Therefore, there must exist a node u £ Br> such that the 
probability that u is visited by time 2T (by at least one of the agents) is less than 1/2. If the adversary locates 
the treasure at u, then the expected time to find the treasure is strictly greater than T, which contradicts the 
assumption. 
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begin 

Each agent performs the following double loop; 
for j from 1 to oo do the stage j defined as follows 
for i from 1 to j do the phase i defined as follows 



Go to a node u chosen uniformly at random among the nodes in Bf, 
Perform a spiral search for time ti := 2 2l+2 /k; 
Return to the source s; 



end 



end 



end 



Algorithm 1: The non-uniform algorithm A^. 




t=2 2i+2 /k 



Figure 1: Illustration of two agents performing phase i. 



at least £ 2 /2 phases i with i > s. Consequently, the probability that none of the k agents finds 
the treasure by time T(€) is at most j~ e I 2 . Hence, the expected running time is at most 

= O (2 s + 2 2s /k) = 0{D + D 2 /k). 
This establishes the theorem. □ 

Fix a constant p>l. We say that the agents have a p-approximation of k, if, initially, each 
agent a receives as input a value k a satisfying k/p < k a < kp. 

Corollary 3.2. Fix a constant p>l. Assume that the agents have a p-approximation of their 
total number. Then, there exists a (non-uniform) search algorithm that is O(l) -competitive. 
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Proof. Each agent a executes Algorithm Ak (see the proof of Theorem 3.1 ) with the parameter k 
equal to k a /p. By the definition of a p-approximation, the only difference between this case and 
the case where the agents know their total number, is that for each agent, the time required to 
perform each spiral search is multiplied by a constant factor of at most p 2 . Therefore, the analysis 
in the proof of Theorem |3.1| remains essentially the same and the running time is increased by a 
multiplicative factor of at most p 2 . □ 



3.2 Unknown Number of Agents 

We now turn our attention to the case of uniform algorithms. 

Theorem 3.3. Consider a function f : N — > N such that ^jLi < 00 • There exists a 

uniform search algorithm that is 0(f(logk))-competitive. 

Proof. Consider the uniform search algorithm -4 U niform described below. Let us analyze the 
performances of the algorithm and show that its expected running time is T(D, k) := <p(k) ■ (D + 
D 2 /k), where 4>{k) :— 0(/(logfc)). We first note that it suffices to prove the statement when 
k < D. Indeed, if k > D, then we may consider only D agents among the k agents and obtain an 
upper bound on the running time of T(D, D), which is less than T(D, k). 



begin 

Each agent performs the following; 
for I from to oo do the big-stage i 
for i from to £ do the stage i 
for j from to i do the phase j 



end 



end 



end 



end 



k 



3 

A, 



2 J ; 



• Go to the node u chosen uniformly at random among the nodes in B{Dij); 

• Perform a spiral search starting at u for time fcy := 2 l+2 //(j); 

• Return to the source; 



Algorithm 2: The uniform algorithm A 



uniform 



Assertion 1. For every integer i, the time until all agents complete big-stage £ is 0{2 l ). 

For the assertion to hold, it is sufficient to prove that stage i in big-stage i takes time 0(2 l ). 
To this end, notice that phase j takes time 0(D itj + 2 l /f{j)) = 0{2 { - l+ ^/ 2 + 2 t /f(j))). Therefore, 
stage i takes time 

O ^(2^)/ 2 + 2V/(. ? ))j =0(2*). 

This establishes Assertion [l] 
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Set s := |~log((D 2 • /(log k)) /k)~\ + 1 and fix an integer i > s. Then, there exists j G {0, . . . , i} 
such that 2 j <k< 2 j+1 . 

Assertion 2. There exists a constant c G (0, 1) such that the probability that none of the agents 
finds the treasure while executing phase j of stage i is at most c. 

To see this, first note that the treasure is inside the ball B(Di.j). Indeed, = J > 

\Jf^j > D. The total number of nodes in B(D itj ) is 0{D 2 i:j ) = 0(2 i+ V/(j)), and at least 
half of the ball of radius y/Uj around the treasure is contained in B(Dij). Consequently, the 
probability for an agent a to choose a node u in a ball of radius y/Uj around the treasure in 
phase j of stage i is 

If this event happens, then the treasure is found during the corresponding spiral search of agent a. 
As a result, there exists a positive constant d such that the probability that none of the k agents 
finds the treasure during phase j of stage i is at most (1 — d ■ 2~ J ) fc < (1 — d ■ 2~ J ) 2J < e~ c . 
This establishes Assertion [2] 

By the time that all agents have completed their respective big-stage s + £, all agents have 
performed fl(£ 2 ) stages i with i > s. By Assertion [2] for each such i, the probability that the 
treasure is not found during stage i is at most c for some constant c < 1. Hence, the probability 
that the treasure is not found during any of those il(£ 2 ) stages is at most l/d e for some constant 
d > 1. Assertion [l] ensures that all agents complete big-stage s + i by time 0(2 s+i ), so the 
expected running time is 0(J2T=o % a+t /d?) = 0(2 S ) = 0(D 2 ■ /(logfc)/fc), as desired. □ 

Setting f(x) := yields the following corollary. 

Corollary 3.4. For every positive constant e, there exists a uniform search algorithm that is 
0(log 1+£ k)- competitive. 

4 Lower Bounds 

4.1 A Tight Lower Bound for Uniform Algorithms 

We say that a function /: N — > N is sub -exponential if / is non-decreasing and there exists a 
constant c < 2 such that f[x + 1) < c • f(x), for every x € N. 

Theorem 4.1. Consider a sub-exponential function f such that 1//0') = 00 • There is no 

uniform search algorithm that is 0(f(logk))-competitive. 

Proof. Suppose that there exists a uniform search algorithm with running time less than f(D, k) := 
{D + D 2 /k)(f)'{k). Hence, if k < D, then f(D,k) < ^!|M j wne re 0(fc) = 2<j>'{k). Suppose, towards 
a contradiction, that <fr(k) = 0(/(logfc)). 

Let T be a (sufficiently large) integer and set D := 2T + 1. That is, for the purpose of the 
proof, we assume that the treasure is actually placed at some faraway distance D greater than 
2T. This means, in particular, that by time 2T the treasure has not been found yet. 

For every integer i < log T/2, set 

ki := ? and A := yjT ■ h/<f>(k z ). 
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Fix an integer i in [l,logT/2], and consider B(Di), the ball of radius Di around the source 
node. We now deal with the case where the algorithm is executed by fcj agents. For every subset 
S of nodes in B(Di), let be the random variable indicating the number of nodes in S that 
were visited by at least one of the fcj agents by time 2T. (We write x( u ) f° r x({ u })-) 

Note that fcj < £>, and, therefore, for each node u <G B(Di) the expected time to visit u is at 
most f(Di, ki) < Df(f>(ki)/ki = T. Thus, by Markov's inequality, the probability that u is visited 
by time 2T is at least 1/2, i.e., Pr(xO) = 1) > 1/2. Hence, E(xO)) > 1/2. 

Now consider an integer i in [2,logT/2] and set Si := B(Di) \ B(Di-i) . By linearity of 
expectation, E(x(5' i )) = ^ MeS E(x(u)) > \Si\/2. Consequently, by time 2T, the expected 
number of nodes in Si that an agent visits is 



n{\Si\/ki) = n 



'A-i(A-A-i) 



= n 



T 



Ah) J ' 



where the second equality follows from the fact that Di — A_i • zcp ^.^ , and the third equality 
follows from the facts that <j>(hi) — 0(f(i)) and / is sub-exponential. 

In other words, for every integer i in [2,logT/2], the expected number of nodes in Si that 

each agent visits by time 2T is £1 ( ^(fe-) ) ■ Since the sets Si are pairwise disjoint, the linearity of 
expectation implies that the expected number of nodes that an agent visits by time 2T is 

'l°g 772 \ /logT/2 



Consequently, X)'°f2 T ^ 2 0(1?) mus t converge as T goes to infinity, and hence so does Y^i=2^ jjjj- 
This contradicts the assumption on /. □ 

Setting f(x) := x yields the following statement. 

Corollary 4.2. There is no uniform search algorithm that is O (log k)- competitive. 



4.2 A Lower Bound for Algorithms using an Approximate Knowledge 
of k 

We now present a lower bound for the competitiveness of search algorithms assuming that agents 
are given approximations of their total number k. As a special case, our lower bound implies 
that for any constant e > 0, if agents are given an estimation k such that k ~ £ < k < k, then the 
competitiveness is fi(log/c). That is, the competitiveness remains logarithmic even for relatively 
good approximations of k. 

Formally, let e : N — > (0, 1]. We say that the agents have a k £ -approximation of k if each agent 
a receives as input an estimation k a for k that satisfies: 

£!-*(£«) < k <k a . 
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(For example, if the agents have a ^-approximation of k, where e is the constant function equal 
to 1/2, then, in particular, this means that if all agents are given the same value k, then the real 
number k of agents satisfies vfc < k < k.) 

Theorem 4.3. Let e: N — > (0,1]. Assume that the agents have a k £ -approximation of k. If 
there exists a <j){k) -competitive algorithm, where <p is a non- decreasing function, then </>(fc) = 
tt{e{k)\ogk). 

Proof. Assume that there is a search algorithm for this case running in time (D + D 2 /k)cf>(k), 
where <fi is non-decreasing. Suppose that all agents receive the same value k, which should serve 
as an estimate for k. Consider an integer W greater than Ak. Set 

T:=2W-m and j := -|— . 

For the purposes of the proof, we assume that the treasure is located at distance D = 2T + 1 , so 
that by time 2T it is guaranteed that no agent finds the treasure. 
For igN, define 

5, := {u : 2 J0+J - 1 < d(u, s) < V 0+l ) . 

Fix an integer i in {J — log k~\ , • • • , log k\ }. Assume for the time being, that the number of 
agents is ki :— 2 2 \ Note that 

k 1 -^ <h< k, 

hence, fcj is a possible candidate for being the real number of agents. Observe that all nodes in 
Si are at distance at most 2- JO+4 from the source, and that \S t \ = e(2 2j0+24 ) = Q(W ■ h). By the 
definition, jo > i + 1. Hence, ki < 2 3o+i ~ 1 < d(u, s), and therefore using the expected running 
time of the algorithm, it follows that for each node u € Si, the expected time until at least one of 
the ki agents covers u is at most 

2d( "' s)2 • m) < ™ ■ m = t. 

Recall that we now consider the case where the algorithm is executed with ki agents. For 
every subset S of nodes of G, let x(S) be the random variable indicating the number of nodes 
in S that were visited by at least one of the ki agents by time 2T. (As before, we write x( u ) 
for %({u}).) By Markov's inequality, the probability that u is visited by at least one of the ki 
agents by time 2T is at least 1/2, i.e., Pr(x(w) = 1) > 1/2. Hence, E(x(«)) > 1/2. By linearity 
of expectation, E(x(Si)) — J2 u es- ^(x( u )) > Consequently, by time 2T, the expected 

number of nodes in Si that a single agent visits is f2(|5j|/fcj) = f2(W). 

Since this holds for any i in { [ 1 ~ 6 2 ^ log k~\ , • • • , |_| log k\ }, and since the sets Si are pairwise 
disjoint, the linearity of expectation implies that the expected number of nodes that a single agent 
visits by time 2T is n(W ■ s(k) log k). Since T = 2W ■ 4>(k), this implies that cf>(k) = fi(e(fe) log k), 
as desired. This concludes the proof of the theorem. □ 



5 Harmonic Search 

The algorithms described in the Section [3] are relatively simple but still require the use of non 
trivial iterations, which may be complex for simple and tiny agents, such as ants. If we relax the 
requirement of bounding the expected running time and demand only that the treasure be found 
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with some low constant probability, then it is possible to avoid one of the loops of the algorithms. 
However, a sequence of iterations still needs to be performed. 

In this section, we propose an extremely simple algorithm, coined the harmonic search 
algorithm^ which does not perform in iterations and is essentially composed of three components: 
(1) choose a random direction and walk in this direction for a distance of d, chosen randomly 
according to a distribution in which the probability of choosing d is roughly inverse proportional 
to d, (2) perform a local search (e.g., a spiral search) for time roughly d 2 , and (3) return to the 
source. It turns out that this extremely simple algorithm has a good probability of quickly finding 
the treasure, if the number of agents is sufficiently large. 

More specifically, the algorithm depends on a positive constant parameter S that is fixed in 
advance and governs the performance of the algorithm. For a node u, let p(u) := d ^ 2+l5 , where c 
is the normalizing factor, defined so that Y^ u eV(G) p( u ) = (Note that c depends on S.) 

begin 

Each agent performs the following three actions; 

1. Go to a node u £ V(G) with probability p(u); 

2. Perform a spiral search for time t(u) := d(u) 2+s ; 

3. Return to the source; 
end 

Algorithm 3: The harmonic search algorithm. 

Using arguments similar to those introduced earlier, e.g., in the proofs of Theorems |3.1| and |3.3| 
one can prove the following result. 

Theorem 5.1. Let 8 G (0,0.8]. For every e > 0, there exists a positive real number a such that 
if the number k of agents is greater than aD s , then with probability at least 1 — e, the expected 
running time of the harmonic algorithm is 0(D + D ^ ). 



6 Conclusion and Discussion 

We first presented an algorithm that assumes that agents have a constant approximation of k and 
runs in optimal expected time, that is, in time 0(D + D 2 /k). We then showed that there exists a 
uniform search algorithm whose competitiveness is slightly more than logarithmic, specifically, 
0(log 1+£ k), for an arbitrarily small positive constant e. We also presented a relatively efficient 
uniform algorithm, namely, the harmonic algorithm, that has and extremely simple structure. 
Our constructions imply that, in the absence of any communication, multiple searchers can still 
potentially perform rather well. On the other hand, our lower bounds imply that to achieve a 
better running time, the searchers must either communicate or utilize some information regarding 
k. In particular, even providing each agent with a fc e -approximation to k (for constant e > 0) 
does not suffice to bring the competitiveness strictly below O(logfc). 

Although the issue of memory is beyond the scope of this paper, our constructions are simple 
and can be implemented using relatively low memory. For example, going in a straight line 
for a distance of d = I 1 can be implemented using O(loglogd) memory bits, by employing a 
randomized counting technique. In addition, our lower bounds provide evidence that in order to 

3 The name harmonic was chosen because of structure resemblances to the celebrated harmonic algorithm for 
the fc-server problem — see, e.g., 1111 . 
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achieve a near-optimal running time, agents must use non-trivial memory size, required merely 
to store the necessary approximation of k. This may be useful to obtain a tradeoff between the 
running time and the memory size of agents. 

From another perspective, it is of course interesting to experimentally verify whether social 
insects engage in search patterns in the plane which resemble the simple uniform algorithms 
specified above, and, in particular, the harmonic algorithm. Two natural candidates are desert 
ants Cataglyphys and honeybees Apis mellifera. First, these species seem to face settings that 
are similar to the one we use. Indeed, they cannot rely on communication during the search 
due to the dispersedness of individuals [3S] and their inability to leave chemical trails (this is 
due to increased pheromone evaporation in the case of the desert ant). Additionally, the task 
of finding the treasure is relevant, as food sources in many cases are indeed relatively rare or 
patchy. Moreover, due to the reasons mentioned in Section [T] finding nearby sources of food 
is of great importance. Second, insects of these species have the behavioral and computational 
capacity to maintain a compass-directed vector flight |14( I36| . measure distance using an internal 
odometer [SSJ , travel to distances taken from a random power law distribution [S3] , and 
engage in spiral or quasi-spiral movement patterns [STJ ISSJ |S3] . These are the main ingredients 
that are needed to perform the algorithms described in this paper. Finally, the search trajectories 
of desert ants have been shown to include two distinguishable sections: a long straight path in a 
given direction emanating from the nest and a second more tortuous path within a small confined 
area |5S[I52]. 
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